Coulombic effect and renormalization in nuclear pairing 
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We investigate effects of the Coulomb force on the nuclear pairing properties by performing 
the Gogny Hartree-Fock-Bogolyubov calculations for the TV = 20, 28, 50, 82 and 126 nuclei. The 
Coulomb force reduces the proton pair energy and the even-odd mass difference by about 25%, except 
for nuclei at and around the proton shell or subshell closure. We then propose a renormalization 
scheme via a reduction factor 7 P for the proton pairing channel. It is found that a single value of 
Jp(= 0.90) well takes account of the Coulombic effect, for nuclei covering wide range of the mass 
number and the neutron excess including the nuclei around the shell or subshell closure. 
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Introduction. The energy density functional (EDF) ap- 
proach (or the mean field approach) provides us with a 
microscopic framework for describing the static and dy- 
namical properties of atomic nuclei from the nucleonic 
degrees of freedom [l| . Because of its numerical feasibil- 
ity, the quasi-local EDF's (i.e. the EDF's represented by 
local densities and currents that include low-order deriva- 
tives, such as the Skyrme EDF) Q have been applied 
to calculations covering wide range of the nuclear chart. 
With the coordinate space representation, the quasi-local 
EDF's are suitable for describing various exotic defor- 
mations 1 and continuum effects of neutron-rich 

nuclei. 

The pairing correlations play a significant role in static 
and dynamic properties of nuclei at low energy HE), it 
is crucially important to construct the pairing channel 
of EDF (pair-EDF) reproducing the pairing properties 
across the nuclear chart 10j. In the approaches using 



the quasi-local EDF's, the pair-EDF is usually taken to 
be local, in order to keep the numerical feasibility. In 
most Skyrme EDF calculations so far, the energy density 
of the form "HP air (r) = A T [p(r)] «*(r) K T (r) (r = p,n) 
has been adopted for the pair-EDF [l|, |ll|, |l2j, where 
P = Pp + Pn is the isoscalar density of nucleons and k t 
the local pair density, by adjusting a few parameters in 
the function A T [p]. It is found that the strength param- 
eter in A p is substantially stronger than that in A n to 



reproduce the observed pairing properties |ll|, [12J . Such 
asymmetry should originate in the dependence of the 
pairing on the neutron excess as well as in the Coulomb 
force which acts only on protons. To include effects 



of the neutron excess, "HP air (r) has been extended as 
•HP air (r) = S T [p(r),pi(r)]«;(r)/s T (r) in Refs. BE}, 
where p\ — p n — p p , though keeping the charge symme- 
try. 

Although the Coulomb force is an important ingredi- 
ent of the nuclear systems, the Coulomb force has not 
explicitly been included in the pair-EDF in most sys- 
tematic calculations because of its non-local nature. It 
was reported that the proton pairing gaps are reduced 
by 20 — 30% if the Coulomb repulsion is treated self- 
consistently 0, E3 ■ It is not likely that the charge sym- 
metric pair-EDF appropriately represents the Coulombic 
effect. An approximate method to take into account the 
Coulombic effect with keeping the local nature could be 
renormalizing the strength parameter of the proton pair- 
EDF as in Ref. 18]. However, it is not obvious whether 
such a simple renormalization scheme works sufficiently 
well. Moreover, an appropriate value of the renormaliza- 
tion parameter and its dependence on Z and N have not 
been known. 

In this paper we investigate Coulombic effect on 
the nuclear pairing by the self-consistent Hartree-Fock- 
Bogolyubov (HFB) calculations, particularly focusing on 
the renormalizability. A numerical method that is ap- 
plicable to wide range of the nuclear chart with a finite- 
range interaction is required for this purpose. Note that 
the HFB theory with a finite-range interaction is practi- 
cally identical to the approach with a non-local EDF. We 
employ the Gaussian expansion method 0, 113] , which 
is adaptable to drip-line nuclei even with finite-range in- 
teractions. For both the particle-hole (ph) and particle- 
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particle (pp) channels, we adopt the Gogny-DlS [2l| plus 
Coulomb interaction with the center-of-mass correction. 
Although we restrict ourselves to the N = 20, 28, 50, 82 
and 126 nuclei, assuming the spherical symmetry, they 
distribute over wide range of the mass number (30 < A < 
220) and the neutron excess (-0.13 < (N-Z)/A < 0.36). 
It is numerically examined whether the Coulombic effect 
on the pairing can be incorporated by a renormalization 
factor for the proton pairing. 

Hamiltonian. We here describe the EDF in terms of 
the effective Hamiltonian. The Hamiltonian for the HFB 
calculations consists of the nuclear part and the Coulomb 
interaction, 



H = H 



N 



(1) 



where H N = K + V N — H cm ' with the kinetic energy 
K, the effective nuclear interaction V N , and the center- 
of-mass Hamiltonian H c rn \ V N may include many-body 
forces, which are often simplified by a density-dependent 
two-body force. The HFB energy can be represented in 
the EDF form owing to Wick's theorem, though includ- 
ing non-local terms in general. For V we adopt the 
D1S parameter-set of the Gogny interaction in this paper, 
which can be employed without the energy cut-off for the 
pairing channel. Since the short-range NN correlation 
hardly influences matrix elements of V c , we use the bare 
Coulomb force for V c . The spherical HFB calculations 
are implemented for the N = 20, 28, 50, 82 and 126 nu- 
clei, by applying the Gaussian expansion method 3, 2^ , 
with the basis functions given in Ref . [i^f ■ It is noted that 
the exchange term of V is treated exactly, and that both 
one- and two-body terms of H c m - are subtracted before 
iteration. 

In correspondence to the expression of the HFB energy 
by the density matrix and the pairing tensor [23| |. we 
separate the Hamiltonian into the pp part H pp that gives 
the pairing tensor and the ph part -ff p h- Each of them is 
comprised of the nuclear and the Coulomb parts; 



H p h — H, 



N 
ph 



V, 



ph ' 



TT _ TTN 



(2) 



We consider the pairing between like nucleons as usual, 
neglecting the proton-neutron pairing, which is not im- 
portant except Z « N cases. H pp is then separable into 
the proton and neutron parts, 



However, the interactions have finite range in general, 
and it is not obvious whether and how the energy of nuclei 
can be approximated to sufficient precision by the local 
limit. In particular, whereas there have been validating 
arguments for iJ p h [2, 24 j, local approximation for H pp 
from which %P alr (r) is derived has not been well explored. 

In the HFB calculations of nuclei, we reasonably pos- 
tulate that H N is isoscalar. Acting only on protons, V c 
breaks the charge symmetry. While the charge symme- 
try is broken at the Hartree-Fock level, there should also 
be difference between H pp for protons and for neutrons 



because of V pp . 



H pp has often been determined so as 



to reproduce the observed pairing properties among neu- 
trons [3, 

To examine whether effects of V pp 
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namely, by using only H pp . 

can be treated in a 
simple renormalization scheme, we define the following 
Hamiltonian, 



H{l P ) = H^+VC+^H^+H^ = (H-Vg)-{l- lp ) H? p , 

(4) 

dropping V pp and introducing the renormalization pa- 
rameter 7 p . While the charge symmetry in the pairing 
channel does not hold because of V pp in H , — (1 — j p ) HP p 
gives the charge symmetry breaking in the pairing chan- 
nel of H(-/p). Note that many HFB calculations so far 
have employed H(l), by presuming the charge symmetry 
for H pp . The central question here is whether or not we 
have 



(H) H « (H( 7p )) R{jp) , 



(5) 



with an appropriate 7 P . We can then renormalize H pp via 
7 p to represent the Coulombic effect. The prescription for 
the Coulombic effect using y v has been applied to 17 Ne 



within a three-body model 27[ ■ Since we here carry out 
self-consistent HFB calculations, the HFB energy at the 
left-hand side (lhs) is evaluated by a calculation with the 
full Hamiltonian H, while the energy at the right-hand 
side (rhs) with H(j p ). These Hamiltonians are explicitly 
written as subscripts in Eq. ([5]). If the wave functions 
are similar, we have (H^ + V p c h + H pp )n ~ (H 
^p P )h(7 p ) an d Eq- © therefore indicates 



V p %)h 



(G) 



which is further reduced to 



pp pp 



H, 



pp ' 



(3) 



(^)h«-((1-7p)^ p }h( 7p)! 



(7) 



The proton pairing should be subject to HP p + V pp . If the 
Hamiltonian contains only the zero-range interactions, 
we need only the local limit of the density matrix and the 
pairing tensor, which leads to a local or quasi-local EDF. 



via {HP p ) H w {HP p ) S ( lp y The value of j p may be de- 
termined for individual nucleus. However, for the renor- 
malization scheme via 7 P to be useful, 7 P has to be fixed 
without referring the result of H for individual nucleus. 
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It is hence desired that 7 P is insensitive to nuclide or ex- 
pressed by a simple function of Z and N. In this work 
we consider the simplest case that 7 P is a constant, with 
no Z or N dependence. 

Results. Figures Q] and [5] depict the spherical HFB re- 
sults of the pair energy (for Z = even nuclei) and the 
even-odd mass difference (for Z = odd nuclei) in the 
N = 20, 28, 50, 82 and 126 isotones 0]. Because the 
neutron pair energy vanishes in these calculations, the 
even-odd mass difference is straightforwardly connected 
to the proton pairing. While several neutron-rich N — 20 
and 28 nuclei seem to be deformed in reality [i^, it is 
sufficiently meaningful to compare the results of H and 
H(lp) within the spherical HFB calculations. The pair 
energy i?P au '(= i<;P air ) is defined by the energy contribu- 
tion of iJpp (i.e. (HP p + V^ p ) H or ( lp HP p ) R ^ p) ), which is 
a simple and clear indicator to the pairing. The even-odd 
mass difference is defined by 

A P (Z) = E(Z, N) - ± [E(Z - 1, N) + E(Z + 1, N)] , (8) 

for Z = odd nuclei. The HFB energies of the Z = odd nu- 
clei are calculated in the equal-filling approximation [3o| | , 

As an observable 
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which has been shown to work well 
directly corresponding to the data, A p has clear physical 
meaning. 



Let us first compare the E^ a " 



and A p values without 



V pp (i.e. by -ff(l), green triangles in Figs. Q] and [2]) and 
those of the full Hamiltonian H (red circles). Similar 
comparison was made in Ref. 17[ , although they viewed 
the pairing gap of the canonical basis locating adjacent 
to the Fermi energy. It is found that, both for £'P air and 
A p , the ratio of the value of H to that of H(l) is about 



75% (Fig. 1 (b,d,f,h,j) and Fig. 2 (b^hj) 
seems consistent with those in Refs. 



dff,hj)) 
IfllllJ. 



. This result 
This ratio is 

almost stable for the nuclides under consideration except 
those in the vicinity of 34 < 42 Si, 90 Zr, 146 > 190 Gd, 208 Pb and 
218 U. 208 Pb is a typical doubly-magic nucleus. The 90 Zr 
and 146 Gd nuclei are well-known as the proton-subshell- 
closed ones, reasonably having suppressed Eg al1 and A p . 
Similar suppression takes place in 34,42 Si, 190 Gd and 218 U 
in the spherical HFB calculation with the D1S interac- 
tion, whereas 42 Si has been suggested to be deformed 
by experiments pjjj . Because of the subshell nature, the 
ground states of these nuclei lie around the boundary be- 
tween the normal fluid and the superfluid phases. Hence 
the usually perturbative V pv affects £"P air and A p to ex- 
ceptional extent. The same consequence was reported in 
Ref. for 90 Zr. 

We next apply the Hamiltonian H(j p ) to the self- 
consistent HFB calculations, adjusting 7p so as to repro- 
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FIG. 1. (Color) Comparison of the pair energy E^ 1 obtained 
from the HFB calculations: (a) E^ lr obtained with H (red 
circles), H(l) (green triangles) and ff(0.90) (blue squares), 
(b) ratios of ££ air obtained with H (red circles) and H(0.90) 
(blue squares) relative to those obtained with H(l), for the 
N = 20 isotones. Analogously, iJ£ alr and ratio for the (c,d) 
N = 28, (e,f) N = 50, (g,h) N = 82, and (i,j) N = 126 
isotones. Lines are drawn to guide eyes. 



H . We find that a single value 7 P = 0.90 satisfies Eq. © 
to good approximation all over the nuclei in this wide 
range of A and (N — Z)/A, as is clear by comparing the 
blue squares with the red circles in Figs. Q] and [2] Re- 
mark that this is true even for the nuclei around 34 ' 42 Si, 
90 Zr, 146 190 Gd, 208 Pb and 218 U, in which the Coulombic 
effect looks exceptionally strong. It has been confirmed 
that the difference in the HFB energies between H and 
fffO.90) is less than 0.1 MeV, indicating that Eq. © it- 
self is fulfilled to good precision. As the wave functions 
of #(0.90) resemble those of H, Eq. Q with 7 p = 0.90 
is good as well. Thus the full Hamiltonian H is well ap- 
proximated by the renormalized Hamiltonian iJ(0.90) in 
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FIG. 2. (Color) Comparison of the even-odd mass difference 
A p and of its ratio obtained from the HFB calculations for 
the (a,b) TV = 20, (c,d) N = 28, (e,f) TV" = 50, (g,h) TV = 
82, and (i,j) TV = 126 isotones. See Fig. [1] for conventions. 
Experimental values taken from Ref. [32J are shown by black 
short-dashed lines in (a,c,e,g,i). 



the HFB calculations, from stable to unstable nuclei in 
wide mass range. 

The purpose of the present work is to investigate the 
Coulombic effect in the pairing channel within the HFB 
framework. However, we also display the experimental 
data of A p for reference in Fig. [5J excluding those of 
the neutron-rich TV = 20 and 28 nuclei that have been 
indicated to be deformed. The HFB results of H and 
-ff (0.90) are in better agreement with the data than those 
of H(l) for the TV = 82 and 126 nuclei, but not for the 
TV = 20 and 28 nuclei. This might imply a room to 
improve the pairing channel in the D1S interaction or 
influence of additional correlations. 

In the present work all the calculations are imple- 
mented by using the Gogny-DlS interaction. However, 



we have confirmed that Eq. (0, in which the lhs is 
the long-range Coulomb interaction while the rhs is a 
short-range nuclear interaction, is well fulfilled for wide 
range of mass region. This indicates that suppression of 
the pair correlation by such weak repulsion is not sensi- 
tive to interaction form. Therefore the renormalization 
of the pairing channel with 7 P will plausibly be appli- 
cable to other interactions or pair-EDF's including lo- 
cal pair-EDF's. Moreover, while the value of 7 P may 
somewhat depend on HgL or the pair-EDF, it should 
not largely deviate from 0.90 as far as the pairing has 
appropriate strength. Although charge-symmetric pair- 
EDF's have been assumed in the usual Skyrme EDF ap- 
proaches the charge symmetry in the pair-EDF 



should be broken because of V£L. It is desirable to read- 
just the pair-EDF (with p and p\ dependence) by taking 
into account this Coulombic effect, for which the renor- 
malization with 7 P (~ 0.90) will be useful. 

Summary. We have investigated influence of the 
Coulomb interaction on the pairing channel in the spher- 
ical HFB calculations. Using the Gogny-DlS plus 
Coulomb interaction for the neutron-closed nuclei of 
TV = 20, 28, 50, 82 and 126, we have found that the 
Coulomb interaction reduces the pair energy and the 
even-odd mass difference by about 25%, except for nuclei 
around the proton-shell- or subshell-closed ones 34 ' 42 Si, 
90 Zr, 146 190 Gd, 208 Pb and 218 U. Because of the non-local 
nature, explicit inclusion of the Coulomb force is not 
adaptable to the local or quasi-local EDF approaches. As 
a renormalization scheme, we have introduced a reduc- 
tion factor for the proton pairing channel of the nuclear 
force (or the pair-EDF), and adjusted the factor to the 
results with the Coulomb interaction. It is found that the 
Coulombic effect is approximated with a single renormal- 
ization factor 7 P (= 0.90) to good precision, all over the 
nuclei under consideration ranging 30 < A < 220 and 
-0.13 < (TV - Z)/A < 0.36, even including the shell- or 
subshell-closed nuclei. 

In the present work we have numerically investigated 
the Coulombic effect on the pairing and the renormaliz- 
ability with H(^ p ). It is of interest to justify the renor- 
malization scheme from a microscopic viewpoint, though 
it is left as a future work. 
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